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1.0  Introduction 

This  report  describes  seme  preliminary  results  of  an  effort  to  relate 
total  aerosol  optical  depth  at  visible  wavelengths  to  the  visibility  in  the 
marine  boundary  layer  (MBL)  .  We  present  an  atmospheric  model  which  explains 
the  observed  height  dependence  of  the  aerosol  extinction,  examine  aircraft 
and  ship  data  frcm  the  MAC^VT  experiment,  illustrate  the  influence  of  synop- 
tic scale  weather  patterns  on  aerosol  structure  and  present  a  simple  model 
to  relate  a     to  total  vertical  depth.  The  ultimate  goal  of  the  project  is 
to  allow  an  estimation  of  IR  ranges  (due  primarily  to  atmospheric  water 
vapor  density,  Q)  from  satellite  remote  sensing  data.  A  feasibility  study 
on  the  concept  was  performed  on  two  NPS  data  sets  (JASIN  and  CEWCCM  78)  us- 
ing the  relationship 

Q/Qs  =  a  +  b  exp  [-1.25  <*o±/2    (7/U)3/2]  (0) 

where  Q  is  the  saturation  water  vapor  density  at  13  C,  <*  the  surface 
extinction  coefficient  in  km"  and  U  is  the  wind  speed  in  m/s.  A  least 
squares  fit  for  834  points  yielded  a  =  1.02,  b  =  -.27  and  a  correlation  of 
0.56. 


2.0  Background 

2.1  Optical  Depth 

Let  the  atmospheric  transmission  coefficient,   T,   due   solely  to 
aerosols  be  represented  by 

T  =  exp  (-u)  (1) 

where 


0O 


(Z)dZ  (2) 

o 
and  a  is  the  aerosol  extinction  coefficient,  Z  the  vertical  coordinate,  and 

/  *r2  dN/dr  C(x,n)dr 
o 
with  r  the  aerosol  particle  radius,  dN/dr  the  aerosol  number  density  spec- 
trum and  Q(x,n)  the  Mie  scattering  efficiency.   If  dN/dr  is  'Known  from  the 
surface  to  the  top  of  the  atmosphere,  then  u  can  be  calculated.  Alter- 
nately, u  can  be  measured  from  satellites.   The  question  is  how  is  u  re- 
lated to  the  surface  visibility  (extinction)  . 

Consider  an  idealization  of  the  atmospheric  extinction  coefficient 
vertical  profile  (Fig.  1) 

u  =  <ct>  n  +  u  +  u  (4a) 

m      t    s 

ut  =  «+  L  (4b) 

where  <<*>  is  the  average  extinction  in  a  well  mixed  boundary  layer  of 
depth,  h,  a  is  the  extinction  in  the  clearer  air  just  above  the  boundary 
layer,  L  is  an  effective  potential  scale  length  for  the  decay  of  a(Z)  for 
Z  >  h,  u  is  the  optical  depth  of  the  troposphere  above  the  mixed  layer  and 
u  is  the  optical  depth  of  the  stratosphere.  This  ideal i2ation  is  a  result 


Extinction 


Figure  1.  Schematic  representation  of  the  height  dependence  of  aerosol 

extinction  coefficient,  a,  within  and  above  the  marine  boundary 
layer  of  depth  h. 


of  the  typical  well-mixed  nature  of  the  marine  boundary  layer   (Section 

2-2).     The   surface  extinction,   a   ,    is  related  to  <a>     by  the  factor  f 

o  m     J 

a     =   f<<*>  (5) 

o  m  v   ' 

which  is  less  than  one  due  to  the  vertical  structure  of  temperature  and 
water  vapor  in  the  marine  boundary  layer. 

If  all  of  the  relevant  variables  can  be  estimated,    then  we  can  calcu- 


late a     using 


aQ  =  f(W  -  \  -  ug)/h  (6) 


For  example,    w     and  w     could  be  obtained  from  the  LOWTRAN  Tropospheric 
Model,  h  frcm  synoptic  climatology  and  f  established  by  model  calculation 
and  empirical  measurement.     Typical  values  of  L  are   2  to  5  kn,    <*    might  be 
10       to  10-1  Im"1  in  the  visible. 
2.2     NPS  Aerosol  Boundary  layer  Model   Structure 

The  marine  boundary  layer  meteorological  structure  is  illustrated  in 
Fig.    2.      The  aerosol   structure  in  the  veil-mixed  context  is  depicted  in 
Fig.    3.      The  total  aerosol  volume  spectral  density,   V(r  ),    is  the  sum  of 
continental,    V  ,   and  sea-salt,   V  ,    components,   where   (V(r  )    =  dV/dr  ) 

V  =  V     +  V  Z  <  h  (7a) 

s         c 

V  =  V  Z  >   h  (7b) 

c 

at  some  reference  humidity  (we  have  assumed  no  clouds)  .  The  volume  density 
is  related  to  the  number  density  by 

V(r)  =  4/3  *r3  dN/dr  (8) 

and  the  particle  radius  at  saturation  S  is  given  by 

r  =  rQ  g(S)  (9) 

The  specification  of  V  and  V  is  based  on  empirical  parameteriza- 
tions.  For  example,  V  can  be  represented  as  a  Junge  distribution  and  V 
as  a  modified  gamma  function,  F  (Goroch  et  al,  1982) 
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Figure  2.     Marine  boundary  layer  extinction  for  virtual  potential  tempera- 
ture,  8   ,  and  water  vapor  density,   Q,   in  the  mixed-layer  model 
structure. 
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Figure  3.     Similar  to  Fig.    2  but  for  the  aerosol  density. 
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Vc  =  A  r^1  (10a) 

Vs  =  VoF(rQ/    S,    y)  (10b) 

where  V  ,    3,  and  Y  depend  en  wind  speed. 

Alternatively/   one  can  use  log-normal  distributions  to  represent  both 
components.      In  a  recent  study  of  atmospheric  electrical  conductivity  this 
method   (which  included  a  log-normal  cloud-droplet  component,   V,)    was  found 
to  be  reasonable .      The  NPS  log-normal  aerosol  number  density  model  is  as 
follows : 

N(r)   =N(r/2^)   exp[-(*n(r/rm)  )2/(2ot2)]  (11a) 

^  =  *n(<0  (lib) 

where  r       =  rm/g(S)    is  the  dry  (S  =  0)   radius 

g(S)  =  0.8   (l-^n(l-S))    ,  S  >  0.4  (12a) 

g(S)   =  1  (12b) 

The  model  parameters  are  given  in  Table  1. 


Table   1 .     Log-normal  aerosol  parameters  for  the  NPS  model .     Wind 

speed  is  u(m/s)   and  cloud  liquid  water  density  is  W(gm/]<g) 
(Note  that  W  depends  on  height  within  cloud.) 


r     ,Mm 
mo 

a 

N 

Cent. 

.027 

2.2 

io2  -  io5 

Sea   Salt 

.088  +  .026  u 

2.2 

0.0012  u3'5/r  J 

'    mo 

CLoud 

4  +  10  W 

1.4 

50     -  200 

A  typical  Junge  value  (A  =  1  at  S  =  .8)    corresponds  to  N     =  650  par- 
ticles per  cm   .     A  typical  wind  speed  of  7  m/sec  gives  55  sea- salt  par- 

3 
tides  per  cm     with  a  number  density  mode  radius  of   0.27  um. 

2.3     Relaxation  of  the  Well-Mixed  Assumption 

Considerable  evidence  exists  to  suggest  that  the  well-mixed  assumption 

is  an  approximation.      It  is  believed  that  this  is  due  to  the  balance  of 
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surface  flux,   entrainment  and  the  boundary  layer  mixing  time.     Ebr  the 

quantity,    X,    let  FQ  be  the  surface  flux  and  F.    be  the  flux  at  the  top  of 

the  boundary  layer  (fluxes  defined  positive  upv\ard)  .     X  can  be  9   ,  q  or 

aerosols.     Define  three  velocities,  W  ,  W    and  W  . 

eg  o 

We  =  "Fh/Ah  (13a) 

WQ=     Fo/Ao  (13b) 

W     =  max(W  ,   W  )  (13c) 

g  o      e 

The  A  quantities  are  the  junps  in  X  at  the  two  interfaces 

4o  =  Xo  -  \o  <14a> 

\  -  Xf  -  Xmh  (14b) 

where  X     is  the  sea-surface  value,   X^  the  value  of  X  at  Z  -   0,    X  ,    the 

value  of  X  within  the  mixed  layer  at  Z  =  h  and  Xf  the  free  atmosphere  value 

of  X  at  Z  =  h  (Fig.   4) .      Ihe  value  of  X  within  the  mixed  layer  at  height  Z 

is  given  by 

X(Z)    =  ^  +  ?(We/2Wg)K+1   Ah  +   (l-C)(Wo/2Wg)K+1   Aq  (15) 

where  5  =  Z/h,    X     is  the  average  "well-mixed"  value  and  K  is  an  integer   (K 

=  3  is  presently  used)  .     NDte  that  W_  =  (L1'2  C  }'2  u  where  C  is  the  drag 

O  2\  w 

-3  .  -3 

coefficient.     A  typical  value   for  W     is  3  x  10       while  W     is  9  x  10       at 

u  =  7  m/s. 

The  application  of  Eqn   (15)   to   9     and  q  leads  to  a  different      (more 
gradual)   height  dependence  of  S  in  the  mixed  layer.      The  application  to 
aerosols  leads  to  a  slight  decrease  in  the  total  number  density  with  in- 
creasing height.      Both  of  these   factors  result  in  a  more  gradual   increase 
of  extinction  with  increasing  height  in  the  mixed  layer.      Therefore,    f  will 
be  closer  to  1.0. 
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Figure  4.      Representation  of  the  gradient  of  a  meteorological  variable 

in  the  boundary  layer  v*hen  the  wll-mixed  assumption  is  relaxed, 
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2.4  Mixed-La ^er  Depth 

The  application  of  Eqn   (6)    is  critically  dependent  upon  a  reliable 
estimate  of  h.      In  raid- latitudes  dominated  by  high  pressure   (e.g.,    U.S. 
West  Coast)   h  is  typically  200  to  800  m,   in  tropical  regions  h  is  typically 
1  to  2.  5  tai.      Note  that  h  can  vary  a  factor  of  two  over  a  24-hour  period  at 
a  given  location.      It  is  hoped  that  h  can  be  determined  from  climatology 
(combined  with  synoptic  scale  classification)   or  remote  sensing  of  nearby 
cloud  top  temperatures. 

2. 5  Synoptic  Classification 

The  method  used  for  synoptic  classification  was  taken  from  the  Refrac- 
tive Effects  Guidebook  (REG,  CGMTHIRFLT  lac.  Mem.  230-2-75).   Figure  5  is  a 
sample  synoptic  chart  frcm  the  REG.   The  REG  classes  are  given  below  in 
Table  2. 


Table  2.   REG  synoptic  classes,  typical  duct  heights,  subsidence 
and  air  mass  type  (ST  =  subtropical,  P  =  polar,  T  = 
tropical,  M  =  modified) . 


TYPE 

h,m 

SUBSIDENCE 

A 

200 

++ 

ST 

B 

600 

+ 

ST 

C 

200 

-H- 

P 

D 

700 

+ 

MP 

E 

1700 

0 

T 

F 

3000 

0 

MT 

G 

4800 

- 

MT 

H 

100, 

1700 

0 

T 

I 

High 

- 

T 

J 

2400 

0 

MT 

K 

no  inversions 

0 

Temp 
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Figure  5.        Sample  synoptic  v*eather  chart  for  the  North  Efecific  in  Spring 
The  capital   letters  designate  regions  of  different  synoptic 
class 
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It  must  be  remembered  that  the  PEG  designations  of  duct  height  are 
based  primarily  on  discontinuities  in  the  water  vapor  profile.     The  lowest 
water  vapor  discontinuity  may  or  may  not  correspond  to  the  top  of  the  mixed 
layer.      For  example,   the  E  category  data   frcm  MAGAT  had  insignificant  water 
vapor  junps  at  the  top  of  the  mixed  layer  (typically  h  -  250  m)  .      The  major 
discontinuity  in  water  vapor  occurred  above  1500  m,   consistent  with  the 
value  given  in  Table   2- 
3-0    MAGAT  Experiment 

The  Marine  Aerosol  Generation  and  Transport  experiment  (MAGAT)  ves 
held  in  Monterey  in  May,    1980.      The  instrunented  platform  participating 
were  the  R/V  ACANIA  (NPS)   and  the  Airborne  Research  Associates   (ARA)   - 
Bellanca  aircraft.      Aerosol   sampling  devices  were   installed  on  both  plat- 
forms ( Eairall ,    1981 ) . 
3.1     REG  Classifications 

Each  day  during  MAGAT   (4/29  -  5/09)  was  classified  using  the  surface 
pressure  chart  analysis  frcm  FNCC.      The  results  are  given  in  Table   3. 

Table  3-      REG  classes  during  MAGAT. 


DATA 

CIASS 

4/29 

E  (Front  2300) 

4/30 

D 

5/01 

E  (Front  2300) 

5/02 

D 

5/03 

D 

5/04 

D 

5/05 

E  (Front  0700) 

5/06 

D 

5/07 

D 

5/08 

E  (Front  1500) 

5/09 

E 
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3. 2  Aerosol  Extinction  Profiles 

Typical  profiles  for  MAGAT  D  and  E  type  classes  are  shown  in  Fig.   6. 
The  MAGAT  profiles  were  divided  into  cloud- free  and  cloudy  profile  subsets. 
A  profile  was  designated  as  cloud- free  if  no  cloud  droplets  ware  encoin- 
tered  during  the  profile  measurements.      This  does  not  imply  that  no  clouds 
ware  visible  anywhere  but  simply  means  that  this  particular  so  aiding  did 
not  pass  through  a  cloud    (Table  4)  . 

The  average  valua  of  f  for  the  cloud- free  cases  was  0.90  with  a  stan- 
dard deviation  of  0.22  and  an  uncertainty  in  the  mean  estimate  of  +  0.06. 
There  was  no  significant  difference  in  the  average  f  for  D  and  E  classes 
but  the  variance   for  D  (based  on  only  four  points)   was  much  less.     Note 
that  the  optical  depth  \i  '    in  Table  4  is  the  total  actually  measured  by  the 
aircraft  (typically  to  a  maximum  altitude  of  1.5  km),   not  the  total  atmo- 
spheric optical  depth.     The  average  mixed  layer  height  under  cloud- free 
conditions  was  h  =  0.29  km  while   under  cloudy  conditions  h  =  0.49  km.     This 
simply  implies  that  the  average  lifting  condensation  level    (LCL)    for  the 
experiment  was  no  greater  than  0.3  km. 

3 . 3  Tropospheric  Extinction 

A  composite  of  all  tropospheric  extinction  profiles  from  MAGAT  (Z  >  h) 
is  given  in  Fig.    7.      Since  the  profiles  do  not  indicate  any  significant 

height  dependence  on  this  scale   (.4  to  1.5  km),    the  data  ware  combined  in  a 

-  -2 

number  distribution  graph  (Fig.   8).      For  D  classes  a     =  1.2  x  10       and   for 

-  —2 

E  classes  a     =  3.0  x  10     .      Tne  standard  deviation  for  LGT   (a.)   was  0.5. 

The  difference  between  a     for  the  two  classes  is  probably  not  due  to  a 

difference  in  aerosol  content  but  a  reflection  of  the  greater  relative 

hunidity  (RH  >  40%)    present  above  the  mixed  layer  for  REG  class  E. 
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Figure  6a.      Perosol  extinction  coefficient   (\  =  0-49  micron)    for  a  D  class 

day 
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Figure  7a.      Tropospheric  extinction   (a  =  0-49  micron)    for  Z   >  h  and  REG 

class  D) 
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Figure  7b.      Tropospheric  extinction   (x  =  0-49  micron)    for  Z   >  h  and  REG 

class  E) 
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Table  4.      Data   summary  for  the   26  ladder  profiles  from  MAGAT. 

p. '    is  the  total  optical  depth  to  the  maximun  altitude 
measured  and  p      is  the  mixed  layer  optical  depth. 

CLEAR  SKY 


Date 

£ 

a     km 
o  — 

yj_ 

-Hm 

h, 

.km 

f_ 

REG 

4/29 

1 

.57 

.33 

.19 

.31 

.92 

E 

4/29 

2 

.46 

.40 

.37 

.55 

.68 

E 

4/30 

4 

.47 

.20 

.19 

.34 

.84 

D 

5/1 

5 

.14 

.074 

.018 

.16 

1.17 

E 

5/1 

6 

.23 

.11 

.082 

.24 

.67 

E 

5/4 

14 

.18 

.065 

.055 

.25 

.82 

D 

5/5 

17 

.17 

.040 

.019 

.15 

1.34 

E 

5/5 

18 

.13 

.035 

.012 

.15 

1.02 

E 

5/5 

20 

.10 

.029 

.020 

.18 

.9 

E 

5/6 

22 

.28 

.18 

.15 

.53 

.98 

D 

5/7 

23 

.46 

.14 

.12 

.23 

.88 

D 

5/9 

26 

.21 

.15 

.088 

.23 

.55 

E 

Average  = 

.28 

.15 
CLOUDY 

.11 

D 

E 

.34 
.25 

.9 

4/30  3  .5  .45  D 

5/2  7  .18  .40  D 

5/2  8  .15  .44  D 

5/2  9  .19  .42  D 

5/3  10  .21  .52  D 

5/3  11  .26  .55  D 

5/3  12  .19  -  D 

5/3  13  .56  .50  D 

5/4  15  .27  .60  D 

5/4  16  .16  .42  D 

5/5  19  .13  .32  E 

5/6  21  .30  .60  D 

5/7  24  .35  .50  D 

5/8  25  .20  .90  E 
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The  estimation  of  the  tropospheric  contribution,  y  ,  requires  a  scale 
length  L  Which  is  either  the  exponential  decay  length  or  the  height  of  the 
tropospheric  aerosol  ccmponent  h  .   If  ws  take  a  typical  valua  of  L  =  7.5 
km,  then  y.  =  0.75  for  class  D  and  p.,  *  .225  for  class  E  for  MAGAT.  Assum- 
ing a  climatological  mean  for  y  =  .01  +  .005  (Toon  and  Pollack,  1976)  then 

s  — 

u   /  y=  .11/. 2  =  .55  for  D  and  p    Ai   =  .11/. 35  =  .32  for  E  cases.      Alterna- 
nt m 

tively,    if  v*e  assume  h    =  2  km  (based  on  REG  data  and  some  aircraft  obser- 
vations)   then  n  /u  =  .11/. 14  =  .79   for  D  and   u  /y  =  .11/ -18  =   .61  for  E 

m  tn 

cases.      Thus,   the  mixed  layer  typically  contributes  60%  to  80%  of  the  total 
optical  depth  for  MAGAT. 
3.4    Variability  Analysis 

The  MAGAT  data  have  yielded  some  estimates  of  the  variability  of  the 
components  in  Eqn   (6)    Which  \*e  can  use   for  a   simple  error  analysis 

(Table  5). 

* 

Table  5.      Variability  contribution  of  the  items  of  Eqn   (6)    for  MAGAT. 


Parameter, 

X 

f 

us 

ut 

h,km 

Average ,    X 

0.9 

0.01 

0.09 

0.3 

AX 

0.2 

0.005 

0.07 

0.15 

Aa     contri . 

o 

Afy  /h 

m 

fAyg/h 

fAut/h 

fy   Ah/h2 

m 

Value   Aa 
o 

0.07 

0.015 

0.21 

0.17 

* 


This  value  of  y,    is  the  estimated  total   y.,    not  the  portion  measured  by 
the  aircraft. 

2 
A  sum  of  the  variances   ( a   )   of  the   four  terms  yields  a  total   standard 

deviation  a  =  0.28  km"     for  a     (therefore  Aa  Va     =  0.85)  . 

3.5     Algorithm  Evaluation 

The  MAGAT  data  can  be  used  to  test  Eqn   (6)    in  a  form  optimized  to  the 

data   set  by  using  MAGAT  mean  parameters   for  the  terms.      If  v\e  assune   a,   = 
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0.01,  H  =  0.35  for  class  D  and  a  =  0.03  and  E  =  0.25  for  class  E,  then 

aQ'  =  .9(u'-  0.013)/1.35  for  class  D  (16a) 

aQ'  =  .9(u'-  0.033)/l.25  for  class  E  (16b) 

In  Table  6  we  show  a  comparison  of  a  '  with  the  measured  a  for  the  12 

Ir^  o  o 

MAGAT  cloud- free  profiles. 


Table  6.      Evaluation  of  Eqn   (6)  as  optimized   (Eqn   (16))    for  the 
mean  conditions  during  MAGAT. 


# 

Class 

u 

m 

y'-u 

aQ' ,   km"1 

a   ,  km 

a7aQ 

1 

E 

.33 

.29 

1.04 

.57 

1.82 

2 

E 

.44 

.40 

1.44 

.46 

3.13 

3 

D 

.20 

.19 

.49 

.47 

1.24 

4 

E 

.074 

•  036 

.13 

.14 

.92 

5 

E 

.11 

.072 

.26 

.23 

1.13 

6 

D 

.265 

.052 

.13 

.18 

.72 

7 

E 

.04 

0 

0 

.17 

0 

8 

E 

.035 

0 

0 

.13 

0 

9 

E 

.029 

0 

0 

.10 

0 

10 

D 

.18 

.17 

.44 

.28 

1.57 

11 

D 

.14 

.13 

.33 

.46 

.72 

12 

E 

.15 

.11 

.40 

.21 

1.9 

An  analysis  of  the  data   in  Table  6  yields  (<*0' /O   =  1-07  and  Aa^a     =   .92, 
which  is  consistent  with  the  estimates  of  the  previous  section. 
4.0     MPS  Model  Simulation 

The  NPS  aerosol  model,  as  described  in  Section  2.2,   has  not  been  run 
on  the  MAGAT  data   sets.     A  preliminary  simulation  was  done   for  a  typical 
climatological  mean  for  open  ocean  conditions.     The   following  parameters 
were  selected:     N    =  650  cm     ,    u  =  7  m/s,   h  =  .35  km,    RH(z  =  10  m)   =  80%, 
Ts  =  15°C,    T  =  13°C,    Aq  =  -5,    A9     =  10  and  W    =  0. 1  an/s.      This  case 
yielded   um  =  0.33  and   f  =  0.72   (Fig.   9)  .      The  value  of  N    was  based  on  open 
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Figure  9.        Sample  NPS  aerosol  model   simulation   for  typical  open  ocean 
climatological  conditions   (see   text,    section   3-4) 
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ocean  measurements  from  CEWCCM-78,    JASIN  and  STOEX.      tfc>te  that  N    =  650 

-3      •  -3 

cm       yields  a     =  3  x  10     ,   which  is  considerably  lower  than  the  values 

found  during  MAGAT.      The  R/v  ACANIA  measurements  of  the  continental  ccm- 

-3 
ponent  imply  a     =  8  x  10       which  is  in  agreement  with  the  aircraft  data. 

The  model  gave  f  =  0.77  which  is  considerably  less  than  the  value  of  0.9 

found  for  the  MAGAT  data . 

5.0  Conclusions 

5.1  Concept 

•  The  optical  depth  is  considered  to  be  made  up  of  three  atmospheric 
regions:  mixed  layer,  tropospheric  and  stratospheric. 

•  The  relationship  of  surface  visibility  (Eqn  6)  is  quite  simple  in 
the  three  layer  context. 

5.2  MAGAT  Evaluation 

The  REG  synoptic  categories  ware  clearly  correlated  with  certain 
features  of  the  extinction  profile   (e.g.    a     and  h)  . 

•  The  value  of  f  of  0.9  indicates  that  it  is  a  smaller  correction 
than  anticipated   frcm  model   simulation. 

The  tropospheric  component  of  the  optical  depth  varied  about  a 
factor  of  three . 

•  Using  parameters  optimized  to  the  data,   Eqn   (6)   could  predict  a     to 
roughly  a  factor  of  2- 

5.3  NPS  Aerosol  Model 

•  The  model  reproduced  the  general  mixed  layer  and  tropospheric 
extinction  profiles . 

Although  the  model  allows  some  relaxation  of  the  well-mixed 
assumption,    it  appears  that  the  mixed  layer  gradients  are  slightly 
greater  than  given  by  the  model . 


The  primary  utility  of  the  model  in  the  future  for  this  project 
may  be   in  its  ability  to  indicate  the  presence  of  a  cloud- free   (h  < 
LCL)   boundary  layer. 
5.4     General   Garments 

•  The  accuracy  of  the  model   (Eqn   (6))    is  primarily  limited  by  the 
variability  of  h  and  u     (see  Table  5) . 

•  Global  application  of  Eqn   (6)   using  climatology  only  will  probably 

allow  estimates  of  a     to  about  a   factor  of   3.      (tote  that  it  is  the 

o 

accuracy  of  the  Wells-Munn-Katz  model  based  on  in-situ  wind  speed 
and  humidity  measurements .) 
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